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We present technical aspects of next-to-next-to-leading order calculations for heavy-to-light decays 
such as top quark decay, semileptonic b quark decay into a u quark, muon decay, and radiative decays 
like b sy. Algebraic reduction of integrals to a set of master integrals is described, methods of 
determining the master integrals are presented, and a complete list of master integrals is given. As 
a sample application, the top quark decay width is calculated to O (os) accuracy. 

PACS numbers: 12.38.Bx,14.65.Ha 


I. INTRODUCTION 

This paper describes the technical aspects of the calculation of O (a^) (next-to-next-to-leading order, or NNLO) 
heavy-to-light quark decays, such as t ^ bW and b ulv, that were presented in [I|. 

Computations in the NNLO for massive charged particle decays are challenging because of the presence of massive 
propagators which complicate loop and phase space diagrams. Here we focus on processes in which we have one 
gluon-radiating particle, Q, in the initial state, and one, q, in the final state, such as a semileptonic decay of a quark, 
Q — !■ 9 -I- leptons. Radiative decays like b —> s"f and b ^ s-\- leptons also fall into this category. An analogous process 
with radiated gluons replaced by photons is the muon decay —> evv. 

A useful approach to such problems consists in expanding the decay amplitudes around some limit in which the 
diagrams can be computed analytically. In the past, all such expansions started with some point on the so-called 
zero-recoil line (see Fig. [T]). Zero-recoil refers to the kinematic case in which the produced charged particle q remains 
at rest in the rest frame of the decaying Q. 

In the present project we address the other kinematic extreme wherein a massive Q decays into two massless 
particles. One of the latter may be a pair of collinear massless leptons or a kF boson in the limit of Q much heavier 
than the W. For a physical top quark decay, we can expand in mw/'mt- Such an expansion, if taken to a sufficiently 
high order, also helps to describe the lepton invariant mass spectrum in the semileptonic b u decay. 

The main technical difficulty in constructing this expansion is the evaluation of a class of diagrams in which the 
four-momenta of all the final-state particles are on the order of the mass of the decaying quark. We explain here 
details of the methods we have used to overcome this obstacle. 

This paper is organized as follows. In Section II we describe the calculational framework: the optical theorem, 
the various contributions needed to construct the expansion, and the ways of computing the necessary integrals. In 
Section III, we list the topologies of the Feynman integrals and give a set of so-called primitive or master integrals to 
which all others can be algebraically reduced. In Section IV, we present explicit examples of how typical primitive 
integrals can be evaluated. Section V contains our results and Section VI summarizes the paper. Throughout the 
presentation, we use top quark decay to illustrate the technical steps. 
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FIG. 1: Kinematic boundaries of the semileptonic decays Q —> g + leptons. The solid arrow shows the expansion presented in 
this paper. Previously known expansions are indicated with dotted arrows Analytical resnlts are now known along the 

whole zero-recoil line [i,i and at the lower left vertex of the triangle (this work). The dot at the left of the triangle indicates 
the kinematic confignration corresponding to the decay t —> hW. 


II. CALCULATIONAL FRAMEWORK 


A. Optical theorem 


At O (q!s), there are three classes of radiative corrections which contribute to the decay t —> bW\ diagrams with 
virtual gluon loops, diagrams with real radiation, and diagrams with both a virtual gluon loop and a real gluon 
emission. The phase space calculations that would be required to calculate the decay rate in this approach are 
difficult even in the case that both masses mb and mw are neglected, and as we mentioned in the introduction, 
finding the functional dependence of the decay rate on mw is one of our objectives. Our task becomes manageable 
by using the optical theorem to relate the decay width to the imaginary part of top quark self-energy diagrams: 


Im(E) 

mt 


( 1 ) 


The self-energy diagrams, E, that we must consider for the O {ctg) decay rate are three-loop diagrams (two gluon 
loops and a t —> bW t loop). The imaginary parts of these diagrams arise from the various cuts that can be made 
through the b, W, and g lines so that the two sides of the cut correspond to particular decay amplitudes. Despite the 
introduction of an additional loop, the calculation is feasible as a result of an assortment of techniques which have 
been developed to solve multiloop integrals analytically with the help of symbolic computation software. 


B. Asymptotic expansions 

With mb = 0, there are two scales in the problem: mt and mw- We define an expansion parameter uj = m'^jmt so 
that the two scales can be expressed as hard and soft {O (1) and O (y/uj), respectively) using mt as the unit of energy. 
The asymptotic expansion in powers and logarithms of to of the contributing integrals is performed using the method 
of asymptotic operation Q , and contributions associated with the two scales are identified using the following heuristic 
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mnemonic. The loop momenta flowing through the gluon lines must be hard or else scaleless integrals (which vanish in 
dimensional regularization) will arise. This leaves us with two regions to consider, depending on the scale of the loop 
momentum flowing through the W line. In the hrst region, all the loop momenta are hard and the Euclidean-space 
W propagator can be expanded as a series of massless propagators: 

1 _ J_ _I I (2) 

k^ k^ 

After the integration, powers of m^/k^ will become powers of uj. 

In the second region, the gluon momenta are hard but the loop momentum flowing through the W is soft. As 
a result, this loop momentum (fc) decouples from the other momenta (£) flowing through the b quark lines via the 
expansion 


I I (P + 2H) (P-P2H)2 {k^ + 2k£f 

[k + £Y~P £4 + Js. +■■■ ■ 

Consequently, the diagrams in this region factor into a product of a two-loop self-energy-type integral and a one-loop 
vacuum bubble integral with a scale of raw- The leading contribution from this second region is O (w^); the interplay 
between the two regions gives rise to terms with a formally large logarithm Inw. 


C. Single-scale integrals 

With the asymptotic expansions in place to extract the functional dependence of the top decay width on mw, we 
have reduced the problem to that of solving single-scale self-energy loop integrals. The indispensable tools for doing 
this are integration-by-parts identities, used to reduce loop integrals with arbitrary powers of denominator factors 
to a small set of master integrals Q. There are two approaches with which these identities can be applied to a 
collection of loop integrals. In the traditional method, one inspects the structure of the identities and rearranges them 
manually into the form of recurrence relations for a complete solution of the system. This method has proven to be 
very successful in numerous applications (e.g., SiE fill, fl^. [T^l 1 but it requires much human work to implement. 

Recently, Laporta [l^ succeeded in developing a targeted solution of the linear system (for only those integrals that 
are needed in a given problem) [Tsj l . This was not previously practical as it is much more exp ensive computationally. 
In our calculation we used the traditional complete approach (programmed in FORM [2) as well as a modified 
version of the new targeted algorithm for which we implemented a dedicated computer algebra program. In both 
cases we independently obtained identical results. This serves as a check of correctness and also enables us to compare 
the two methods El- 


III. LIST OF TOPOLOGIES AND MASTER INTEGRALS 

In the “soft” loop-momentum region, the diagrams factor into a product of a two-loop self-energy-type integrals 
and a one-loop vacuum bubble integral with a scale of mw- The bubble integral is purely real and can be solved 
easily. The entire set of two-loop self-energy integrals was first solved in [T^; the contributions from this region are 
therefore quite straightforward. 

In the “hard” loop-momentum region, we need to extract the imaginary parts of genuine three-loop integrals. We 
identify 9 master topologies in terms of which all 36 O (a^) diagrams can be expressed. The diagrams which give 
rise to these master topologies are shown in Fig. [2J They are labeled by the letters A through I in analogy with the 
notation of El for the O (a^) contributions to muon decay. We use dimensional regularization (H = 4 — 2e) and take 
the top quark mass as the unit of energy, mt = I. Using the notation \d^k] = d^k/{2TT)^ we define the following 
Euclidean-space loop integrals {p is the incoming momentum). 


A(ai, 02,03, 04, 05,06, 07, 08, ag) = 


5(01,02,03,04,05,06,07,08,09) = 


[d^ ki][d^ k2][d^ ks] 


ki (fci -I- p)^“2 (fc^ -1- A:3 -I- p) 2“3 (fcj _|_ ^2 -I- A:3 -I- p)2“4 

I 

^ [(fc2 + ks+p)^ + I]“Kfci + 2A:3p)“6fc2a7^2as(2fep)“9 

[d^k,][d^k2][d^k3] 


I kl^^^ki 


-I- p^2a2 -j- k2 + p)^“3 (^ki +k2+ ks + p)2“4 


( 4 ) 
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FIG. 2: The nine C?(Qs) diagrams for t —> bW which constitute the master three-loop topologies of the calculation. 

( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 

(9) 


<7(01,02,03,04,05,06,07,08,09) = 


£>(01,02,03,04,05,06,07,08,09) = 


£(01,02,03,04,05,06,07,08,09) = 


£(01,02,03,04,05,06,07,08,09) = 


<7(01,02,03,04,05,06,07,08,09) = 


[(^2 + ks+p)'^ + l]“5(fc| -h 2A:3p)“6fc2“"fc3“»(2A:2p)“9 

_ [d^ ki][d^ k2][d^ ks] _ 

kl^^iki + p)2“2(fc3 +k2+ p)2“3fc2“'‘(fc2 + 2k2p)‘^^ 

1 

[(^2 + ks +p)'^ + l]“®(fc| + 2k^pY'! kl°‘^ {2kik^Y^ 
_ [d^ ki][d^ k2\[d^ kz] _ 

kY^ (fci -I- p) 2“2 (^3 -h + p)^“^ (fci + ^2 + ^3 + t)^“'‘ 

1 

^ (fci -h + pY'^^kY^kYYkl + 2k3pYY‘^k2pY^ 

_ [d^ ki][d^ k2][d^ ks] _ 

kYYki +pY'^Hki + k2+pY°-Yki + fca -hp)2“4fc2“5 

1 

^ k^Yk^ + 2k2pYYkl + 2k3pYY‘^k2k3Y^ 

_ [d^ki][d^k2][d^k3] _ 

kY^ (fci -I- ^2 + ^3 + p)^“^ (fci + fe + p)^“3 (fci -I- A:3 -I- p)2“4 

1 

^ kY^kYYkl + 2k2pYYkl + 2k3pYY‘^k2k3Y^ 

_ M^fei][d^fc2][rf^fe3] _ 

^ 1 “^ (fci +k2+ pY°-‘^ (fci -I- + A:3 -I- p)2“3 (/;2 _|_ 2k2pY^ 
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[(fc 2 + ks+p)^ + l]“5(fc2 + 2/c3]9)“6fc2“"fc^“«(2A:i}9)“9 

f [d^ ki][d^ k2][d^ ks] 

-H (01,02,03, 04, as, 06,07, 08,09) = / ,2a, m ^ -i- , 2 a„, 2 a ./,9 -^ 

J /cf l(fci + k2 +p)2“2(fc^ + /C3 +p)2“3/c2a4;,2a5^^2 _^2fc2p)“6 


(fc| + 2fc3p)“^(fc2 — A:3)2“8(2 A:ip)“9 


/(ai, 02 , 03 , 04 , 05 , 06 , 07 , 03 , 09 ) = j ^ 2 ai^^ 


[ii^fci][(i^fc 2 ][(i^fc 3 ] 


+ p)2“2(fc^ -(- k2 +p)2“3(fc| + 2k2p)°‘^{kl + 2fc3p)“3 

1 


kl°‘^kl'^^{k2 - fc3)2“«(2A:iA:3)“® 


( 10 ) 


( 11 ) 


( 12 ) 


For the decay rate calculation, we need the imaginary parts of these integrals, obtained by setting p^ = — 1 + zO. 
For the exponent 09 we only need non-positive values; it refers to a product of momenta that may remain in the 
numerator after canceling scalar products against denominator factors. 

In order to specify the boundary conditions of the recurrence relations, we need 24 master integrals. They vary 
greatly in their complexity, but for completeness, we shall list them all. Some of these integrals — particularly those 
with six or fewer lines — are used by more than one of the topologies and thus any descriptions we provide using the 
integral labels A through / are not necessarily unique. JF denotes the loop factor r(l -|- e)/(47r)'°/^. 

Five of the master integrals can be solved in closed form to all orders in e via a sequence of simple one-loop integrals: 


Im 


[d^k 


[d^k 


_ _ [ Kks] 

kliki + pY j (fc| + 1)7 (fc| + 1) 


Im 


f [d^ks] 

kfk^iki + k2 + pY J (fc| + 1 ) 


[d^ki][d^k2 


= 


= 


I 4 

“2 ■*" “ 

c 


47r"^ 

II-— ) + ( 26---2^3 ) e 


+ (57-i|^-8C3 + ^)c^ + 0 (e^) 


1 15 

^ ~ T ' 

1155 157r2 


145 TT^ 


+ SCa 1 e 


Im 


[d^ki][d^k2Wk^\ 
klk^k^ki + k2 + ks + pY 


= 


+ - 


I 


16 4 

8281 , 1457r2 ^ 75C3 7tt'^ 

'W ^ 8 ^ 

71 


(13) 


(14) 


0{e^) 


12 72 

■( 


/3115 TT^ 


432 


f 109403 7l7r2 7(3 




Im 


[d^ki][d^k2][d^k3] 


Im 


[d^k 


(fci + pY{ki + k2Y{ki + ksYk^kl 

f [d^fc2][d''A:3] 


kl{ki+pYJ fc|fc|[(fci + ^ 2 +p )2 + 1 ] 


= 


= 


2592 36 

/ 3386467 3II57r2 497(3 477^ 

V 15552 216 1^ I5 

+ 10 + (64 - 277^) e -p (336 - 2077^ - 22(3) 

77r4 


(15) 


+ 1584 - 12877^ - 220(3 + 


6 


1 


2e2 4e 


47 TT^ 


“^72 “ 77 + “"S-^ + —77- A -3(3 e 


o{Y) 


133 37r2 


(16) 


16 


V417_W_^_^\ , 

V 32 24 2 4 y '' '' 


(17) 


Two more master integrals are also a sequence of one-loop integrals, except that the final one-loop integral is less 
trivial: 


Im 


[d^k3][d^k2][d^k3] 


[{ki +pY + l](fci + k2Yiki + *3 +pYklk'l 


= 


7 7r2\ /175 37r2 

7 - TT + -^ - IIC 3 e 


(18) 


/2681 13l7r2 99(3 2977^ 

24 2 ^ 


+ O {Y) 
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Im 


[d^ ki][d^ k2][d^ k^] 


k\[{ki +p)^ + l](fci + fc2 + fca +p)2A:|fc| 


= 


1 25 


383 




(19) 


One of the master integrals requires the real part of the two-loop massive self-energy integral first evaluated in 


Im 






kl{ki -\- pY J (fc| -I- l)(fc| -I- 1)[(A:2 + k3+ pY -I- 1] 


= 


3 29 

765 137r2 


175 


TT 

~2 


+ - 


1943 977r2 


12 

27C3 tt 


3C3 


( 20 ) 


32 


24 


2 -^ + 8.^1n2)e^+0(e3) 


Five of the master integrals are only needed to leading order in our calculations and are very closely related to some 
of the four-loop integrals calculated for muon decay in M- Specifically, if the massless neutrino-electron loop is 
integrated out, an overall factor of 5(1 — e, 1 — e)iF/e is obtained and the exponent on the W line increases from 1 
to (1 -|- e). To leading order in e, neither this modification to the W exponent nor the difference in the normalization 
conventions changes the master integral, and so after converting the appropriate results from Minkowski space to 
Euclidean space, we conveniently obtain the following master integrals: 


Im 71(1,1,0,1,1,1,1,1,0) 
Im 5(1,1,1,1,1,1,1,1,-1) 
Im 5(1, 0,1,1,1,1,1,1,0) 
Im5(l,l,l,l,0,0,l,l,0) 
Im 5(1,1,1,1,1,1,1,1,-1) 


7r5^ 




177r4 


+ Oie) 




13<3 


' 4 


-io 



+ 0{e) , 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


The other integrals in [l^ are also relevant to our calculation, however, since we require analytic results beyond the 
leading order in e for the remaining master integrals, we can only use these four-loop integrals as a consistency check 
on our own calculations. The remaining eleven master integrals are: 


Im 


Im 


Im 


[d^ki][d^k2\[d^k3\ 


[kl -I- l](A:i -I- fc2)^fc|fc|(fc2 + k3+ pY 


[d"fci][d"fc2][d"fc3] 


kf{ki -I- pYk'^k'^Kki + k2 + k3 + pY + 1] 


[d^ ki][d^ k2][d^ k3] 


(fci +pY{ki + k2Y{ki + k3Ykl{kl + 1) 


= 7r5^ 


= 7r5^ 


= 7r5^ 


- 5C3 e' + O (e") 


1 15 /145 TT^ 

^ T [~8 2 

(1155 IStt^ 

1 5 17 / 49 , 

+ (-™ + 10<3) «2 + 0 (,») 

M /II 7r2\ (11 37r2 ^ \ 

/439 197r2 ^ 


(26) 


(27) 


■ (28) 


Im 5(1,1,1,0,1,1,0,1,0) 
Im 5(1,1,1,1,1,2,0,0,0) 


TT^'^ 


TT^'^ 


^ + ^+ (?---2C3 




131 77r2 77r4 , ^^, 3 . 


+ ^^ -I- TT^ In2 -I- — In^ 2 - - In'^ 2 - 8Li4(l/2)^ e + O (e^) 


, (29) 
(30) 
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Im 5(1,1,1,1, 2, 2, 0,0,0) 

Im (7(1,0,1,1,0,1,1,1,0) 
Im (7(1,1,1, 0,1,1,1, 0,0) 
Im F(1,1,0,0,1,1,1,1,0) 

Im F(l, 1,0,1,0,1,1,1,0) 
Im (7(1,1,1, 0,1,1,1, 0,0) 
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3C3 , 


TT' 

"6 4 ' 2 

TT^ 3(^3 IItt^ 
_ 


In 2 


+ TT^ In 2 + y In" 2 - i In^ 2 - 8Li4(l/2) ) e" + (7 (e^) 


1 + (14 - 4C3) e + 119 - tt" - 28C3 - 




1 5 

-1 + ^ 


15 TT^ , 
y-Y-2C3) + 


Ttt^ 


29 tt" 
~2~~Q~ 


0{e^) 


+ -14 + — +9C3 e 


6 


+ -119- 


IStt" 

2 

^2 


63C3 


4l7r4 


ll-y-C3 


77- 


137r" 


6 


+ O (e^) 


- 13C3 - 




+ (7(e") 


TT 


C3 


14-^-3C3 + ^)e + (7(e") 


(31) 

5 

(32) 

(33) 

(34) 

(35) 

(36) 


The integrals in Eqs. (l30l) and (1311) contain the fourth-order Polylogarithm Li4(l/2) = 0.517479..along with terms 
with higher powers of In 2. None of these constants appear in our final results, which suggests that there might exist a 
basis of master integrals in which these constants do not arise in any of the terms at the orders in e which contribute 
to our final results. 


IV. EVALUATION OF MASTER INTEGRALS 


The eleven master integrals in Eqs. (1^ through (15^ were evaluated on a case-by-case basis by integrating the 
phase space of massless propagators severed by cuts. In this section, we will illustrate the procedure with a few 
examples. We will compute master integrals having two-, three-, and four-particle cuts, respectively. 

The general procedure is as follows. For any particular integral, we can use the optical theorem to obtain its 
imaginary part via 

Im A4 = i ^ / dHeut MiM *2 , (37) 

cuts 


where A4i and A42 are the scalar amplitudes arising from the uncut propagators on the two sides of the cut, the cut 
propagators are assumed to be on mass shell, and f dllcut is the integral over the phase space of the cut propagators. 
Working in dimensional regularization, the phase space integral involves a factor of 


(27r)^-i25 


(38) 


for every particle involved in a given cut, along with a 5-dimensional i5-function to conserve overall four-momentum. 
Much of the phase space integration can be done for a general case without assuming anything about the amplitudes 
A4i and A42, so that the phase space integrals relevant to an n-particle cut can be reduced to simpler integrals 
involving a small number of scalar parameters. We shall see this explicitly in the following three examples, each with 
a different number of cut particle lines. 


A. Example 1: two-particle cut 


Our first example will be the integral of Eq. dsni). As can be seen in the sketch of Fig. [31 there is a single cut 
(through the two rightmost propagators) that can be made through massless lines in the diagram. The two-body 
phase space integral for massless particles is especially simple, since the two particles must emerge back-to-back while 
partitioning the total energy equally. Defining the general phase-space factor by 


n 


jD-1, 


DxD 


Li (27^) 


D-1 


2K; 


(27r)^6 



(39) 
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FIG. 3: Sketch of the master integral 13(1,1,1,1,1,2,0,0,0). The imaginary part of this integral can be obtained from a 
two-particle cut through the massless lines on the right side of the diagram. Solid and dashed lines depict massive and massless 
propagators, respectively. 


an explicit calculation leads to the two-body phase space factor 


T 22'=7r3/2 

r(l-b e)r(3/2 - e) ’ 


(40) 


where, as before, T denotes the loop factor r(l -b e)/(47r)^/^. There are no remaining integrals in this expression, 
and so we have 

ImB(l,l,l,l,l,2,0,0,0) = ^ Li , (41) 

where Li is the factor arising from the two-loop integral for the three-point function to the left of the cut in Fig. [S] 
Explicitly, 


Li 


_ [d°k][d^l] _ 

{k +pY[{k +p - qY + 1][(A: -b lY + 


(42) 


where k and I are the loop momenta, p is the incoming momentum from the left (with the onshell condition p^ = — 1), 
q emerges onto the massless line on the upper-right (g^ = 0), and (p — q) emerges onto the massless line on the 
lower-right (such that {p — qY = 0). Applying the Feynman parameters x, y, and z in three consecutive stages, we 
can evaluate the I- and /c-integrals so that we are left with the parametric integral 


Li = 


r(l-b2e) 
r2(l + e) 





xYl- x)-^z{l- z f 
[z'^xy + zx{l — 2y) -b 1 — 


(43) 


The y-integral can be solved explicitly, but the remaining x- and 2 -integrals are more difficult. The integral is finite 
(as a result of our use of a second factor of one of the massive lines) and thus we can expand the integrand in e using 
expressions like 


= 1 -b e In a; -b ^ In^ a; -b O (e^) 


(44) 


and then evaluate the integrals order by order in e. The leading (e°) term is relatively simple, but the presence 
of products of logarithms in the subsequent terms quickly escalates the level of difficulty. The FORM packages 
Summer and Harmpol have proven to be very useful for evaluating some of these integrals; some of the 
related formulas are conveniently tabulated in [2^. For the first two terms of Li, we obtain 


Li=T 


2 




SCa IItT^ n n 1 4 / , X 

-b--b7r2ln2-b — ln^2- - In^ 2 - SLU 1/2) | e- 

2 240 3 3 / ;I 




(45) 


and, in conjunction with Eqs. (1401) and (1411) . thereby establish the result in Eq. dsni). 


B. Example 2: three-particle cut 

An example of a master integral requiring a three-particle cut is Eq. (1551) . This integral is sketched in Fig. 21 where 
it is seen that a diagonal cut through the massless lines leaves a one-loop amplitude to the left of the cut and a 
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FIG. 4: Sketch of the master integral 1^(1,1, 0,1, 0,1,1,1, 0). The imaginary part of this integral can be obtained from a 
three-particle cut through the massless lines from the upper left side of the diagram towards the lower right. 


single massive propagator to the right. The three massless particles that we cut through can be put on shell with 
a range of relative kinematic configurations, meaning that the three-body phase space integral can only be partially 
solved in advance. The remainder of the integral will be over two kinematic parameters (such as the energies of two 
of the particles) upon which the amplitudes M.\ and M .2 might also depend. A very convenient formulation of the 
three-body phase space involves the use of the relativistic invariants 

u = {li + (l-iY and 2 =(<72+ < 73 )^ (46) 

as the integration variables. This leads to the expression 

fP '\2 /*! pl — z 

P 3 = ^ p {I - z- u)~^ , (47) 

Z-kBii Jo Jo 

where Bn = 73(1 —e, 1 —e). Returning to the master integral in Fig.[4l we observe that the massive propagator (m = 1) 
in the upper right part of the diagram splits directly into two onshell massless particles, and thus the momentum q 
flowing through this propagator is related to one of our invariants: 

(g2 + l)-i^(l-z)-i . (48) 


In addition, we must include a factor in the phase space integral resulting from the one-loop integral in the lower left 
portion of Fig. |4l The momentum I flowing through this diagram is related to the other invariant via P = —u, and 
the loop integral can quickly be evaluated with a Feynman parameter: 

/ [(t +‘P+ w ° i 

Substituting this expression into Eq. dUD, we have 


Im 7^(1,1,0,1,0,1,1,1,0) 


AttcBii 


(50) 


where J denotes the 3-parameter integral: 


J = 



(1 - x)-’^ u-^ z-^ {l-z- u)-^ (1 - xu)-^ 


(51) 


As in the previous example, this integral is finite, and hence we can safely expand the integrand in e and evaluate it 
term by term. The result. 


J — 1 -|- 



e-b 



IIC 3 



e^ + 0{e^) , 


(52) 


once folded in with the other factors in Eq. (EOl), leads to the expression of Eq. (ESI). 


C. Example 3: four-particle cut 


Our final example illustrates how a four-particle cut can be used to obtain Eq. (IM)) . This integral is sketched in 
Fig. [SI On each side of the cut we have a massive propagator. As might be expected when cutting through four lines 
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FIG. 5: Sketch of the master integral 1^(1,1, 0,0,1,1,1,1, 0). The imaginary part of this integral can be obtained from a 
fonr-particle cut through the massless lines. 


in a planar three-loop diagram, there are no loop integrals remaining. The most general expression for a four-body 
phase space integral is extremely complicated, but it happens that we can obtain a simpler expression suitable for our 
purposes by modifying the three-body result. Specifically, the generalization of the three-body phase space integral in 
Eq. (H71) from the case of three massless particles to the case of two massless particles and a particle of squared-mass 
y = q 2 (such that 0 < y < 1) leads to 


Psm = 


(P2)^ 

2ttBii 


/•i 

dz du u 

Jo Jo 


{1- z){z-y) 



(53) 


When y = 0, this reduces to Eq. dUD, but if we interpret y as the relativistic invariant of a pair of massless particles 
in a four-body phase space. 


2 / = (<?2 + qi )"^ , 


(54) 


we obtain a four-body phase space integral 
(p,)3 ri ri- 


/'z{l — z — u)/{l — z) 


Jo Jo 


du 


dy u 


y-^ 


[I- z){z-y) 


(55) 


that is applicable to certain integrals where two of the massless particles form a subloop. In order to evaluate the 
integral in Fig. 0 we will need to replace the invariant u = (gi -|- qzY with v = {qi + 92 )^ since the massless loop 
corresponds to a massive q 2 and thus the massive propagators beside the cut carry momenta {qi + q 2 ) and (92 + 93 )- 


With 


u+v+z=l+y, 


(56) 


Eq. (1551) can be rewritten as 


P 4 = 


(P2)^ 

‘iir'^B 


r-i+y-z 


11 Jo 


dy / dz 


'vh 


dv y ^z ^{l + y — z — vY 



(57) 


The integral we need to evaluate has the additional propagator factors (1 — 'c) ^ and (1 — z) and as before, can be 
expanded in e and evaluated term by term, thereby leading to the result of Eq. (l34l) . 


V. RESULTS 


The QCD corrections to the decay width for t bW can be written as 


r(t ^ bW) = To 


Xo + ^X. + i^)" X2 + 0{al) 

TT V TT / 


where 


To = 


Gf \Vtb\^ wY, 


8\/27r 

is the tree-level decay width in the limit that mw = 0. The results for Xq and Xi are known analytically 


(58) 

(59) 
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The two-loop contribution with which we are presently concerned, X2, has been estimated numerically 
Detailed comparisons of our analytical results with those studies is given in [l[ . We divide X 2 into four gauge-invariant 
color structures, 


X2 = Cf {TrNlXl + TrNhXh + CfXa + CaXna) ■ 


(60) 


Here we use SU{N) factors (with N = 3) Cp = {N‘^ — l)/(2iV), Ca = N, and Tr = 1/2. Np denotes the number 
of massless quark species and Np that of quarks with the mass equal to that of the decaying quark. Dehning the 
parameter to = our results for the four classes of contributions in Eq. ((6OI) are: 


Xl = 




+UJ° 


+LU‘ 


237r2 

108 

5839 

689 



19 27r2' 


-b OJ 

7~^ ^ 

-b 


27000 
2282381 
'37044000 


Ca 

' 

TT^ 

15 ~ 900 


2C3 

7 


\nuj 


TT 

70 


3 

Inw 

2263 

88200 


+ UJ^ 


+ iO° 


Inta 


745 

4253 
8^ ' 
13187 
181440 ’ 


SlTT^ 7 ■ 


TT 

T ’ 

TT^ 

' ^ 


17 

144 

1 

' 48 


Inw 


In w 


(61) 


Xh = 


12991 537r^ 


1296 


54 


3 


35 

108 


47r^ 


dCa 


6377 257r^ 


432 


18 


+iO'= 


■319 

3l7r^ 

2Cal 

-bw'* 

'76873 

00 

to 


'237107 

87r21 

. 27 

27 

3 

8640 

9 

27000 

9 


+ Ca 

+ O (tu^) 


(62) 


= 


5- 


IIOtt^ 53C3 


+UJ 




-bw^ 


-ba;° 


48 

7537 


8 720 

5237r2 295C3 


19 2, 

—TT^ In 2 
4 


19l7r4 


73 4l7r2 

-TT + 


288 

16499 


96 

40771^ 


32 720 

7Ca , Ttt^ 


-b LU 

27 2, 

-TT^ In 2 - 

16 


115 


4177^' 

90 

57r^ 


864 216 

1586479 295171^ 


^„ , 367 57r2 , 

2 ' 120“' >"2+1-143 + ^''"“ 


9 


259200 
11808733 
' 6480000 


6912 

377r2 


2400 5 


2 

6Ca 


-b 


-In a; 

17280 16 


/13589 TT^ , , 

-Inw 

I 27000 60 


In w 


+ 0{u;^) , 


(63) 


Xna = 


521 


+UJ^ 


+uj-= 


+UJ^ 


+UJ^ 


5057r^ 

864 

12169 

576 


16 


llTT^ 19 2, „ 

-TT^ In 2 

1440 8 


217l7r2 377C3 


7771-* 


- UJ 


27 


91 , 3297r2 
48 144 


137r' 


41 


576 


64 


288 32 


+ TTTrTT ln2 -b 


73 

16 


60 

37r^ 


In a; 


13685 477r2 I9C3 4377^ 1 


864 72 4 

420749 326371^ 


-b 


720 2 

9C3 /11941 


-b -7r^ln2-b - 


1121 


432 


TT 

"6 


In a; 


37r^ 


103680 13824 

4868261 5577r2 


6912 32 


Int 


12960000 4800 


10 


/153397 TT^ , , 
i 216000 ~ 40 


+ 0(w®) . 


(64) 


As discussed in [l|, we obtain the two-loop QCD correction to the top quark decay width by substituting w ~ 0.213 
into these expressions. There are already sufficiently many terms so as to render the theoretical uncertainty 20 times 
smaller than the experimental uncertainty induced by the mass measurements of the top quark. Furthermore, these 
expressions can be smoothly matched to the corresponding expressions Q for the two-loop QCD corrections to the 
semileptonic decay b uiv in the zero recoil limit. The result of such a matching procedure is depicted in Fig. [SI 
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X2 



FIG. 6: Matching of expansions around a; = 0 (thick line) and oj = 1 (thin line). The solid line denotes the resulting correction 
to the decay width valid in the fnll range of ui. 


VI. CONCLUSIONS 

We have presented an analytical approach to next-to-next-to-leading order studies of heavy-to-light decays. We 
have used top quark decay as an illustrative example here; in [ij we also studied the distribution of the lepton 
invariant mass in the semileptonic 6 — > it decay and have reproduced the O (a^) corrections to the muon lifetime. A 
number of other phenomenologically interesting studies become possible with the results presented here. For example, 
corrections to radiative b sj decays and to Bd,s meson mixing require matrix elements of very similar if not identical 
kinematic structure. (For the mixing, one can expand in the four-momentum carried by the light quark and obtain 
propagator-like diagrams studied in the present paper.) 

We employed two approaches to reducing Feynman integrals to a basis set of primitive ones. In both we use 
integration-by-parts recurrence relations that we either solve completely, or for a targeted set of needed integrals. 

The most appealing feature of the targeted approach is that the core of the program is universal for all topologies. 
Unlike the complete solution method, where each new topology necessitates a by-hand construction of a new solution 
algorithm, the targeted approach needs only a minimal specification of the topology under consideration. As a result, 
once the program core is written and tested, very little programming is needed to incorporate new topologies, thereby 
saving time, effort, and minimizing the likelihood of programming bugs. 

Another major difference between the two approaches is the use of subtopologies. In the traditional method, 
the removal of a line will often facilitate a mapping of a loop integral from a more complex topology onto a simpler 
topology, for which a new set of identities is used. The drawback to this intuitively reasonable process is a proliferation 
in the number of topologies to be programmed. For example, the 9 master topologies that we encounter in this work 
(discussed in section IIIII) require approximately 40 additional subtopologies. The targeted approach, on the other 
hand, only needs to deal with the master topologies, as it is able to solve a given topology completely using only one 
set of identities. At the three-loop level, it would be very difficult for a human to program an algorithm in this way. 

A key feature of the targeted approach is the generation of large tables containing all possible integrals in a topology 
with a predetermined range of powers. Once this table is generated, it is accessed by the programs responsible for 
assembling a particular physical calculation. In the complete method, conversely, the algorithms are typically applied 
in real-time, and the generation of dynamic tables will likely be used in the near future to improve the performance 
of this approach. Given the large range of powers which can arise when performing a calculation using asymptotic 
expansions, these tables and intermediate expressions can become quite large. 

In summary, the targeted approach appears to be superior by virtue of saving human time. The burden is shifted 
to the computer and thus an efficient implementation of the algorithm is essential. We developed the project-specific 
software based on the BEAR package [^. The targeted approach has also recently been implemented in MAPLE 

M- 

The availability of the targeted approach greatly simplifies multiloop calculations. The remaining bottleneck involves 
the evaluation of master integrals. Eurther progress in perturbative calculations may depend on a radically new 
approach addressing this issue. 
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APPENDIX 


As described in Section El expansion around w = 0 is given by the functions dSHl-il. The complementary 
expansion around lo = 1 has been studied in [J. Below, the corresponding formulae are given for completeness: 


+(5^ 

+(5^ 

+(5® 

Xna{S ~ 0 ) = 

+(5^ 

+5" 


6C3 + 


117 41 

~+18^ 

■797 53 ^ . 

108 27^ C3+(^ 18 

'1289 1 , 37, . 1 


39 4 . 

31 8 


ln(5- -In^J 
2 

ln(5- iln^5 

3 


432 18 

133 5 2 

~ 6^ 


—TT^-In (5 H— In S 


18 


5400 30^ 15 


6 

797 1 2 

3^ 

+ O (<5®) 


■ 

r1817 3 J 

+ (5® 

-In J 


[3600 5 J 




1 


5400 36 10 


(A.l) 


(A.2) 


523 71 2 7 2 , „ 84 39, 

-TT^ + -TT^ In 2 H-TT^-C 3 - 

32 12 2 15 4 ^ 


1363 217 2 1 2 , o 32 4 


147 

le 


— 27r^ ) ln5 H-In^ 6 




lnJ+ i^ln2 J 

V 24 

9 y 

4 


1537 

937 2 

5 2,,. 

/67 

17 2 ^ 



-TT^ - 

-TT^ In 2 — 

-f 

— TT^ 

In J -I-fn^ d 

288 ' 

864 

12 

V12 

36 y 

8 J 


43609 7609 2 1 . 

43200 21600^ ~ 4^^ ~ 10' 


Tn2 


/269 

8 2^ 

. 11 2 

1 — + 


In (5 + — In^ (5 

V 90 

90 y 

6 J 


, (A.3) 


[llOO 

881 2 

7 2 , „ 

13 4 

129, / 

423 

47 2 ^ 

l , J. 33 2 r 


- TT^ - 

—TT^ ln 2 + 

— TT^ - 


- u 

— TT^ 

1 In (5 H - In S 

[ 32 

144 

4 

60 

8 ^ V 

16 

12 y 

8 J 


623 247 


-1" 

54 

- TT^ + 

48 

18319 

19 2 

1728 

572^ 

52379 

5041 

43200 

43200' 


1 2 , „ 13 4 

-TT^ In 2 -TT^ 

6 90 


I 2 C 3 + ( ^-27r" )ln5+^ln"(5 


24 


^ 1^2. 


/1877 

1 2 ^ 


1- 

-TT^ 

in d -In^ d 

V 288 

24 y 

24 J 


40 


C3 


/139 

1 2 ^ 

. 11 2 / 

1-f 

-TT^ 

In (5 - In^ 6 

V 90 

90 y 

48 


0{S^) ,(A.4) 


where 6 = 1 — lo. These two sets of expressions can be smoothly matched in the intermediate region to produce 
functions valid in the full range of u. For practical purposes it is convenient to replace this complicated combination 
of exact expansions by a collection of polynomial fits with numerical coefficients. We found that the following simple 
functions 


Xl 2.854 - 0.665 w-0.109 a;2_g 572 a;3 +5.561 a;^ +0.931 , (A.5) 

Xh -0.063615 +0.098146a;+ 0.144642^2- 0.307331a;3 + 0.107417w^ +0.020707(A.6) 
A ,4 ~ 3.575- 2.867a; + 2.241a;2- 12.027a;3 + 11.564a;^ - 2.489a;^ , (A.7) 

AAr ,4 ~ -8.151 + 2.990a;-3.537a;2 + 36.561 a;3- 42.275a;^ + 23.899a;®-9.494a;® , (A.8) 

give a very good approximation. For any value of to they deviate from the exact value by less than 0.01. The second 
order contribution to the decay rate for = 5 and Nh = 1 can be written in the approximate form 

A 2 ~ -16.729+ 3.064 a; + 4.481 a;2 +35.191 a;® - 10.336 a;'‘- 15.655 a;®, (A.9) 


where the error introduced by the fit is smaller than the uncertainty of the result due to, for example, higher-order 
QCD corrections. 
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